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Abstract. We describe the effects of both magnetic buoy¬ 
ancy and differential rotation on a disc of isothermal gas 
embedded in a purely azimuthal magnetic field, in order 
to study the evolution and interplay of Parker and shear¬ 
ing instabilities. 

We perform a linear analysis of the evolution of perturba¬ 
tions in the shearing sheet model. Both instabilities occur 
on the slow MHD branch of the dispersion relation, and 
can affect the same waves. We put a stress on the natural 
polarization properties of the slow MHD waves to get a 
better understanding of the physics involved. The mecha¬ 
nism of the shearing instability is described in details. 
Differential rotation can transiently stabilize slow MHD 
waves with a vertical wavelength longer than the scale 
height of the disc, against the Parker instability. Waves 
with a vertical wavelength shorter than the scale height of 
the disc are subject to both the Parker and the transient 
shearing instabilities. They occur in different ranges of ra¬ 
dial wavenumbers, i.e. at different times in the shearing 
evolution ; these ranges can overlap or, on the contrary, 
be separated by a phase of wave-like oscillations, depend¬ 
ing on the strength of differential rotation. 

These analytical results, obtained in a WKB approxima¬ 
tion, are found to be in excellent agreement with numerical 
solutions of the full set of linearized equations. Our results 
can be applied to both galactic and accretion discs. 

Key words: accretion discs - instabilities - MHD - cos¬ 
mic rays - Galaxies: ISM - Galaxies: magnetic fields 


1. Introduction 

The Parker instability (Parker 1966) relies on the natu¬ 
ral tendency of the horizontal component of the magnetic 
flux to escape towards the surface of a vertically stratified 
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atmosphere (Kruskal & Schwarzschild 1954, Tserkovnikov 
1960, Newcomb 1961, Yu 1966). This effect of magnetic 
buoyancy occurs wherever a magnetic field is embedded 
in a density gradient: its classical manifestations are the 
sunspots on the surface of the sun (Parker 1955), the es¬ 
cape of the magnetic field lines towards the halo of disc 
galaxies (Parker 1966, Lachieze-Rey et al. 1980), or the 
formation of the molecular clouds which appear as “beads 
on a string” in galaxies (Parker 1966, Mouschovias et al. 
1974, Blitz & Shu 1980, Elmegreen 1982b). 

The latter application, however, has been criticized 
(Zweibel & Kulsrud 1975, Parker 1975b, Balbus 1988), 
and recent observations indicate that the Jeans instabil¬ 
ity is a major ingredient in this process (Elmegreen et al. 
1994). It may be faster than the Parker instability, and 
self-gravity has the advantage of not relying on the lit¬ 
tle known topology of the magnetic field lines. The self 
gravitational collapse may even be helped by the presence 
of a sheared magnetic field (Lynden-Bell 1966, Elmegreen 
& Elmegreen 1983, Elmegreen 1987). Some purely grav¬ 
itational mechanisms have also been advocated (Balbus 
1988), as well as thermal instabilities (Elmegreen 1989a,b, 
1990). Nevertheless, the Parker mechanism may still play 
a role (Elmegreen 1991, Hanawa et al. 1992), especially 
since recent work by Giz & Shu (1993) shows that, with 
more realistic vertical equilibrium profiles, one should ex¬ 
pect a growth rate a factor two higher than in classical 
calculations. The Parker instability could also be a key 
ingredient of the dynamo process in stars (Parker 1975a) 
and galaxies (Ruzmaikin et al. 1988, Hanasz & Lesch 1993 
and references therein), if combined with differential rota¬ 
tion. 

The recent re-discovery of a powerful shearing insta¬ 
bility (Velikhov 1959, Ghandrasekhar 1960) in magnetized 
discs by Balbus & Hawley (1991) has renewed the inter¬ 
est taken in differentially rotating magnetized discs: an 
axisymmetric instability was found to occur in any mag¬ 
netic configuration which is not purely azimuthal (Balbus 
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& Hawley 1991, Dubrulle & Knobloch 1993, Kumar et al. 
1994). It could be efficient even in weakly ionized discs 
(Blaes & Balbus 1994). A vertical magnetic field is also 
unstable to non-axisymmetric perturbations (Tagger et al. 
1992, 1993, 1995). 

The configuration with a purely azimuthal magnetic 
field is stable to axisymmetric perturbations, but tran¬ 
siently unstable to non-axisymmetric ones. This different 
instability was very briefly mentioned by Acheson (1978) 
in the context of stellar interiors, and described in more 
detail by Foglizzo & Tagger (1991), and Balbus & Hawley 
(1992b) when applied to galaxies and accretion discs. 

One of the most promising results of these shearing 
instabilities, combined with the Parker instability, could 
be the description, in a single model, of both the turbu¬ 
lent viscosity and the dynamo action in accretion discs 
(Vishniac & Diamond 1992, Tout & Pringle 1992). Unfor¬ 
tunately, these attempts are still mainly qualitative, be¬ 
cause of the analytical difficulties introduced by differen¬ 
tial rotation. 

On the other hand, although the case of an initially 
vertical magnetic field leads to the most powerful shear¬ 
ing instability, its nonlinear development gives birth to a 
strong azimuthal magnetic component (Zhang et al. 1994). 
This azimuthal magnetic held is subject to both the Parker 
and the transient shearing instabilities, and has not been 
fully studied yet: Balbus & Hawley (1992b) canceled the 
effects of magnetic buoyancy by studying the shearing in¬ 
stability in the disc midplane. Moreover, they introduced 
a radial component of the initial magnetic field, thus pre¬ 
cluding the possibility of an initial hydrostatic equilib¬ 
rium. 

In this paper we describe in detail the effect of differen¬ 
tial rotation on waves in a disc embedded in an azimuthal 
magnetic field, at the linear stage. This has already been 
the purpose of an earlier study (Foglizzo & Tagger 1994, 
hereafter [FT]), where we considered perturbations with 
a vanishing vertical wavenumber because this is known to 
be the most favourable configuration to the Parker insta¬ 
bility, though it considerably limits the effect of shear. 
Here, by considering perturbations with a non vanishing 
vertical wavenumber, acting on a true hydrostatic equilib¬ 
rium (with a purely azimuthal initial magnetic held), in 
a vertical density gradient, we can look at the details of 
the conjugate effects of both the Parker and the azimuthal 
shearing instability. We will show that they occur on the 
same branch (the “slow wave”) of the MHD dispersion re¬ 
lation, at different values of the radial wavenumber. Since 
shearing motions lead the latter to vary with time, this 
means that the same wave can evolve through periods of 
stable oscillations, and shearing or Parker amplification, 
corresponding to a changing polarization (i.e. the orienta¬ 


tion of the perturbed velocities associated with the wave), 
respectively more radial or more vertical. These two am¬ 
plification mechanisms are found to act successively or 
sometimes simultaneously, depending on the strength of 
differential rotation. 

Our paper is organized as follows: in Sect. 2 we first 
recall our hypotheses, which essentially reduce to the 
“canonical” ones used in classical works on the Parker 
instability and to the use of the “shearing sheet” model 
used in the description of spiral waves in differentially ro¬ 
tating discs. This leads to our set of linearized differential 
equations. In Sect. 3 we discuss in detail the validity of 
the WKB approximation, and interpret it in terms of the 
natural polarization of the slow MHD waves. Throughout 
the paper we will And that, in addition to giving a clear 
understanding of the physics involved, the WKB approx¬ 
imation gives results in excellent agreement with exact 
numerical solutions. 

In Sect. 4 we explain the mechanism of the shearing in¬ 
stability. We compare the characteristic features of the 
shearing and the Parker instabilities, taken separately. We 
show analytically in Sect. 5 how these two instabilities 
proceed from the same slow MHD mode, and solve the 
apparent contradiction with the transient shearing stabi¬ 
lization found in [FT]. We also derive the overall strength 
of the transient amplification. 

The phase of transition between the two instabilities is 
shown to depend crucially on the strength of differential 
rotation in Sect. 6 . We summarize our main conclusions 
in Sect. 7, and give the details of our analytical treatment 
in the appendices. 

2. The equations 

2.1. Hydrostatic equilibrium 

The equilibrium is exactly the same as in [FT], and we re¬ 
call it here very briefly. In the rotating frame at distance 
To from the rotation axis, we define the vectors x along 
the radial direction (x = r — rg), y along the azimuthal 
magnetic field (Bq = Boy) and z parallel to the rotation 
axis (Jl = flz). 

The vertical gravitational field is taken to be constant 
(—sign( 2 ;)( 7 zz). This classical hypothesis allows an easier 
mathematical formulation. 

The magnetic pressure (and the cosmic ray pressure in the 
case of a galactic disc) are assumed to be proportional to 
the thermal gas pressure: 

Bb = aPth and Per = /37th- (1) 

The gas is convectively neutral (7 = 1), and both the 
sound speed a and the Alfven speed Va = a(2a)^/^ are 
independent of height. 
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Thus the vertical density profile of the hydrostatic equi¬ 
librium is simply: 

Po{z) = po(0) exp(-|z|/i7), (2) 

where the scale height is: 

= (1 + a+ /?) — . (3) 

9z 

In our numerical calculations, we assume that the vertical 
gravity and rotation frequency are related by: 

9z^—9r^H9?. (4) 

ro 

This neglects the gravity of the midplane disc of stars, and 
consequently would overestimate the scale height H in the 
case of a realistic galactic disc. 

2.2. Linearized perturbations in the shearing-sheet approx¬ 
imation 

We will write our set of equations in terms of the La- 
grangian displacement vector dehned by: 

(^ + V<,.V)^ = v + (eV)V„. (5) 

The perturbations of the density and magnetic field take 
the simple form: 

p = -V.ipoO, (6) 

b = V X (e X B„). (7) 

In the shearing sheet approximation, differential rotation 
is measured by the Oort constant A: 


Defining ^{kx,t) by Eq. (10) permits us to interpret one 
effect of differential rotation as a linear time-dependence 
of the radial wavenumber of the sheared perturbation 
{Kx{t) = kx — 2Akyt). This appears on the differential 
system satisfied by ^{kx,t)- 

0 = X f - (11) 

where the dimensionless hermitian matrix C{kx) includes 
additional terms due to the vertical wavenumber kz' 

^ kl-k2a{kl-kkl-k\) kxky {^-i-P-ikz{l-i-2a))ikx^ 

^ {-^ + Oi+P-ikz)iky 

^ X X 2cKA:y-t-(A:^-j- — ) (1-t-2o() j 

We have noted k = Hk the dimensionless wavenumbers 
normalized by the scale height of the disc, and crosses 
denote the complex conjugate of the symmetric terms. 
The differential force appears through the dimensionless 
parameter kA defined as: 




AAflH^ 


> 0 . 


( 12 ) 


As stated in [FT], the shearing time-scale Tshear = 
occurs both: 

(i) as the typical time-scale of the shearing of waves, 
i.e. the linear growth of the radial wavenumber: 


Kx{t) = kx - 2Akyt, 


(13) 


ro 9D 

.4 = -^(.„) < 0. 


( 8 ) 


In the frame rotating at radius ro with the rotation fre¬ 
quency D = f^(?'o)j tbe sheared motion is linearized so that 
the azimuthal speed of the gas at equilibrium is: 


(ii) as a scaling of the radial differential force, which is 
opposite to the radial magnetic tension. Their ratio is: 


^ 2n Tlarker 
k^Vl kl ^ kl ToTshear 


(14) 


Vo = roDo -I- 2Ax. (9) 

The only difference with [FT] is the presence here of a 
vertical wavenumber kz , so that perturbed quantities vary 
as: 

^{x,y,z,t) = 

_ 1 ^+oo 

p^^e^Vyy+k.z) / 

^ — oo 

kz was set to zero in [FT] for the sake of simplicity, 
since classically this is the most unstable Parker mode. In 
Eq. (10) and hereafter the (fc^, A:^)-dependence of ^{kx,t) 
is not written explicitly, for the sake of clarity, since linear 
calculations allow us to keep them constant. 


3. WKB approximation 

3.1. Mathematical transformations 

The presence of a vertical wavenumber kz still allows the 
mathematical transformations performed in [FT]: we re¬ 
duce the sixth order differential system & to a second 
order differential equation, by first performing a Laplace 
transform: 

^{kx,u;)=e (15) 

In Appendix A we derive an ordinary differential equa¬ 
tion, of which f(kx,oj) is the unique solution bounded at 
infinity. This equation has a singularity where A(cu^) = 0, 
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and this polynomial is interpreted in [FT] as the asymp¬ 
totic dispersion relation at \Kx{t)\ —> oo. 

Reciprocally, is determined by the inverse Laplace 

transform: 


1 i-ip+co . , \ _ 

J 20—00 


(16) 


f ip—OQ 


where p is real and larger than the largest imaginary part 

of the singularities of _ 

The solution of the differential equation (A14) can be writ¬ 
ten in the WKB approximation: 


1 Ap+ca 

iikx,t) = — (/i+e*'^+-Ai_e*'^-)da;. (17) 

Jip—oo 

and the functions p.±{ijj) vary slowly if a WKB criterion, 
derived in Appendix B, is satisfied. We show that the ma¬ 
jor contributions to the integral (17) are then given by the 
saddle points uj±{t) defined by: 

^(u;±(t)) = 0. (18) 

As in [FT], we can move the integration contour of 
Eq. (17) so as to pass through these saddle points, and 
follow the steepest descent path between them. When the 
WKB criterion is satisfied, this formulation allows us to 
interpret Eq. ( p^ as a dispersion relation whose roots are 
time-dependent. The saddle points are the relevant instan¬ 
taneous frequencies or growth rates which used to be the 
poles of the Laplace transform ^(w) in the absence of dif¬ 
ferential rotation. 


3.2. WKB validity and the linear coupling of waves 

It was shown in [ET] that if the WKB criterion is satis¬ 
fied, the evolution is “adiabatic”, i.e. a wave (stable or not) 
preserves its “identity” throughout its sheared evolution. 
On the other hand, at low \Kx{t)\, the WKB approxima¬ 
tion fails if the differential rotation is realistic, and a lin¬ 
ear coupling of waves occurs: an initial pure Parker mode 
(slow MHD wave) is associated with a complex pattern 
of motions, due to the interplay of the Coriolis force and 
magnetic tension with its basic perturbed velocities. As 
a result it loses its identity and it can emit both magne- 
tosonic (fast MHD) and Alfvenic waves in the disc. 

The WKB criterion can be understood as the require¬ 
ment that the instantaneous frequency and polarization of 
the waves evolve (because of the time dependence of the 
radial wavenumber Kx{t)) sufficiently slowly with time. 
This criterion is satisfied in two different cases (see Ap¬ 
pendix B): 

(i) If Kx{t) varies very slowly with time, the very slow 
temporal variations of the growth rate and polarization 
allow a globally adiabatic evolution of the solution. This 
is made possible if the shearing parameter A is very small. 


The effect of the differential force {AAVl^x) becomes negli¬ 
gible compared to the radial magnetic tension {kyV^f^x) if 
the magnetic field is not very weak (as in Shu 1974). Nev¬ 
ertheless, the case of a vanishing field allows us to study 
the effect of the differential force on the Parker instability 
([ET] and here, Sect. 5). 

(ii) If the polarization and the frequency of the pertur¬ 
bation depend weakly on the radial wavenumber Kx(t), 
the WKB analysis is valid even if the shearing parameter 
is not small. We study in the next section the polarization 
and the frequency of the slow MHD wave to determine the 
range of wavenumbers Kx , k^ that allows this approxima¬ 
tion. 

3.3. The polarization of the slow MHD mode 

The slow MHD mode is characterized by displacements 
with a strong component along the magnetic field lines, 
taking advantage of the anisotropy of the magnetic pres¬ 
sure. The amounts of radial and vertical displacements 
depend crucially on the orientation of the wavevector. In¬ 
dependently of the presence of any instability^ the vertical 
stratification and the rotation of the magnetized plasma 
are responsible for the following properties: 

^ = O ^ oo when \kx\ > k^,Hky 

^ = O f ^ ^ 0 when \kz\ > kx, k’^af^/Tl 

\^Z J 

^ ~ when |fc^|, Ifcx] > 

The transition between these asymptotic polarizations oc¬ 
curs at kx ^ kz (if kz > Hky, kyo'^/D,). It should be re¬ 
membered that in most cases of interest for us, Hky ~ 
kyjn ~ i/H 

These properties are typical of both the Alfven and slow 
MHD modes in the classical case of a stratified plasma 
with uniform rotation, and can be derived in a straight¬ 
forward way by a normal mode analysis (see Appendix C). 
The WKB approximation, which is valid in each of these 
asymptotic limits, allows us to generalize these properties 
to the case of differential rotation. 

We know from [FT] that differential rotation causes 
a wave with an initial radial wavenumber Kx ^ — oo to 
evolve linearly with time towards Kx -l-oo. A slow MHD 
wave with an initial Kx{t = 0) <C — [fcz | is polarized mostly 
in the vertical plane {y, z) (see Eq. (19)). It rocks naturally 
towards the horizontal (x, y) plane when, because of the 
shear, \Kx{t)\ < kz (see Eq. (20)). It rocks again towards 
the vertical {y, z) plane when Kx{t) > kz (Eq. (19) again). 
The intrinsic nature of the slow MHD wave is consequently 
enough to explain that the displacements are successively 
vertical, horizontal, and vertical again, and thus more 
prone to buoyancy or Coriolis forces. Figure Q illustrates 


(19) 

( 20 ) 
( 21 ) 




T. Foglizzo et al.: The Parker-Shearing instability in azimuthaly magnetized discs 


5 


qualitatively this evolution. 


Fig. 1. Polarization of the slow magnetosonic mode during its 
sheared evolution. 


When the radial wavenumber \Kx\ ^ kz,Hky, the po¬ 
larization of the slow MHD wave is essentially vertical and 
the eigenfrequency remains constant. The dispersion rela¬ 
tion reduces to A{uj^) = 0, which is independent of kx- 
The evolution is then adiabatic and allows a WKB study 
of the asymptotic Parker instability (Shu 1974, [FT]). 

When the vertical wavenumber \kz\ ^ kx, kyal/^, 
the polarization is essentially horizontal and the eigen¬ 
frequency remains constant. According to Eq. (|^), the 
duration of the phase of horizontal polarization scales as: 

Thoriz. ~ 4 ^- (22) 

Aky 

Thus the polarization of the slow MHD wave is con¬ 
stantly horizontal on an arbitrary long time if the ver¬ 
tical wavenumber is high enough. The dispersion relation 
reduces to Q^(u;^) = 0, which is independent of kx (see 
Appendix B). The evolution is then adiabatic and allows 
a WKB approximation. 

The strongest effect of the differential force is expected 
in this phase of horizontal motions, since this force is radial 
and proportional to the radial displacement. Indeed, if the 
differential force dominates the radial magnetic tension, 
the shearing instability will result, during a time compa¬ 
rable to Thoriz. (see Sect. 4-5). 

In the opposite case where fcj < kyol/^ ~ 1/77, the 
phase of horizontal polarization shrinks to about one shear 


time, and the evolution may be adiabatic only in the limit 
of vanishing differential rotation (see [FT] for kz = 0). 


If both the vertical and radial wavenumbers are large 
Kx ^ kz ^ Hky,kyal/fl, the polarization of the slow 
MHD wave varies slowly: 


dKx{t)/kz 

dt 


-2A^ < |A|, 

kz 


(23) 


and the WKB approximation is possible again (see 
Sect. 6). 


4. The pure shearing instability 

4 ..I. The pure shearing instability 

If a magnetic field line of length Ay = 2'K/ky is bent 
by a radial displacement ^x, a magnetic tension {Fmt = 
PokyV^^x) acts as a restoring force. The differential force 
{Fd = AAVlpQ^x) is also proportional to the radial displace¬ 
ment, but in the opposite direction. Hence if the following 
criterion is satisfied: 

4AfI -b klVl < 0, (24) 

a radially bent field line will naturally continue its radial 
displacement. 

This balance is indeed a central feature of the shearing 
instability. However, this simple description is not enough 
to fully explain it: it ignores the effects of the other forces, 
namely the Coriolis force, the pressure force and the other 
components of the magnetic force. In particular, the Cori¬ 
olis force is usually stabilizing by transforming any radial 
motion into an epicyclic oscillation, as in the case of sound 
waves modified by rotation. 

Moreover, one would like to understand why the criterion 
( P^ seems to predict an instability in the limit of a vanish¬ 
ing magnetic field, which seems to contradict the stability 
ensured by the Rayleigh criterion (k^ = 4D^ -b 4Ail > 0) 
for a non-magnetized disc. 

We now proceed to a detailed description of the mech¬ 
anism, by analysing the linearized system dO) in the “hor¬ 
izontal regime” where kz 1/77, kx, and Eq. (20) is valid. 
In particular the stratification and cosmic ray effects dis¬ 
appear here. 

An illustrative picture of the mechanism is given in Fig. 

It is convenient, for the sake of clarity, to decompose 
the magnetic force into two contributions, which have sim¬ 
ple expressions in our approximation: 

(i) an anisotropic “pressure” force acting against the 
compression of the field lines. 


f Kx 

Fpb = ipoVlV±(y4j_) = -ipoVl{Kx^x + kz^z) d 0 (25) 

\kz 
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Fig. 2. The shearing instability viewed from above the disc: the 
rotation of the disc is here anti-clockwise. The amplitude of the 
perturbation is exaggerated for the sake of clarity. The circles 
along the perturbed magnetic field line represent the maxima 
of density. When the displacements of the slow MHD mode are 
mainly horizontal (Az < Ai,), the radial thermal pressure can¬ 
cels the anisotropic magnetic pressure. Instability occurs if the 
radial magnetic tension is dominated by the differential force. 
The azimuthal thermal pressure force acts as a restoring force, 
which can overcome the stabilizing influence of the Coriolis 
force. 


(ii) a magnetic tension which can be visualized as the 
result of the “elasticity” of the field lines when they are 
bent. 

Ux 

Ftb = -ipoVlklU = -^PoVlkl \ 0 (26) 

The z—component of Eq. (0) implies that in a slow 
MHD mode in the horizontal regime, the thermal pressure 
force cancels the magnetic anisotropic pressure force, to 
first order: 

V(?V.(a) + a'V.^^0 (27) 

Note that this does not mean incompressibility: Eq. ( p^ 
can be equivalently written in order to express the density 
perturbation as a function of the azimuthal displacement: 


P_ _ ikyVl ^ 
Po a^s+Vr^ 


( 28 ) 


Independently of the presence of any instability, this im¬ 
plies in particular that the extrema of the density per¬ 
turbed by a slow MHD mode occur where the azimuthal 
displacement vanishes (see Eig. ||). 


(i) As a consequence of Eq. (|^), the sum of the ther¬ 
mal pressure force and the magnetic force, in the radial 
direction, is reduced to the single magnetic tension term 
mentioned above. 

(ii) As a consequence of Eq. (|^), the azimuthal com¬ 
ponent of the thermal pressure force can be expressed as 
a restoring force acting on the azimuthal displacement. 

Hence in the limit of small vertical wavelength, the 
slow MHD mode (the Alfven mode as well) corresponds 
to perturbations satisfying the following 2-D system: 

The only remaining forces are the Coriolis force, the ra¬ 
dial magnetic tension modified by the differential force, 
and an azimuthal restoring force (we have noted 14 = 
asVx/{a1 -\- the cusp speed). 

We can notice the similarity with the case of the axisym- 
metric instability of a disc embedded in a magnetic field 
which is not purely azimuthal: Balbus & Hawley (1992a) 
obtained an equation similar to Eq. (29), with the follow¬ 
ing differences: 

(i) kz appears instead of ky here, because of the differ¬ 
ent geometry of the equilibrium magnetic field, 

(ii) the azimuthal restoring force involves the Alfven 
speed instead of the cusp speed here. 

They stressed the link with the Hill equations, and gave a 
simple mechanical analogue: the motion of a particle at¬ 
tached with a spring to a guiding center in orbit in a grav¬ 
itational potential. The spring is isotropic (C^ = Cy = 
klV^)- This analogy also holds in our case, by considering 
an anisotropic spring {Cx = kyV^,Cy = kyV^). One can 
be easily convinced that the motion of the particle is un¬ 
stable if the stabilizing effect of the spring is canceled by 
the differential force (4AH -|- Ca, < 0 ). 

The importance of the azimuthal restoring force in the in¬ 
stability mechanism is now clear: the Coriolis force would 
be able to reduce the shearing instability to a neutral mode 
if the azimuthal spring constant Cy were to vanish (for ex¬ 
ample if the sound speed Os —> 0 ). 

This occurrence of the cusp speed may be traced back to 
the fact that we are here dealing with the destabilization 
of the slow MHD mode, whereas the axisymmetric insta¬ 
bility studied by Balbus and Hawley corresponds to the 
destabilization of the Alfven mode. 

As shown in Appendix B, the WKB criterion is sat¬ 
isfied in the horizontal regime, so that we can replace 
time derivatives in Eq. (29) by simple (—*w) multiplica¬ 
tions (note that radial derivatives have disappeared). The 
resulting dispersion relation reduces to Qa^"^) = 0 : 

4-[4+ 2{l+a)klV^]4 + klV^AAn + klVl) = 0 (30) 
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According to this dispersion relation, instability indeed 
occurs if the differential force dominates the magnetic ten¬ 
sion. 

The other root of this dispersion relation, corresponding 
to the Alfven mode, is always stable and reduces to an 
epicyclic oscillation (inertial wave) in the case of a vanish¬ 
ing magnetic field. 

The missing mode is of course the fast magnetosonic mode, 
corresponding to fast vertical oscillations, which was ex¬ 
cluded from our approximation by Eq. (20). 

As the Parker instability, the shearing instability re¬ 
quires a long enough wavelength along the magnetic field 
lines: in the galactic case, the minimum azimuthal wave¬ 
length scales as t:Va\AVI\~^/'^ ~ l.Skpc. 

The optimal azimuthal wavenumber corresponds to 
kivl = (1 + a)ujg — 2Afi, where Ws is the optimal growth 
rate: 


® 1 ■ 

(1 -f 2 q:) 2 -|- 2q; -|- 2q:q) ^ -f 1 -|- 2 q; -f qQ: 

jwsl is a decreasing function of the intensity of the mag¬ 
netic field a. Thus the highest growth rate corresponds to 
a vanishing magnetic field, and we obtain in this limit, for 
an infinite optimal kyi 

^ -A^. (32) 

This extreme limit is of course never reached, because the 
diffusion of the magnetic field must be taken into account 
as soon as the magnetic field is weak: our idealized MHD 
equations then break down, as shown in details by Ache- 
son (1978). 

Nevertheless, this calculation confirms the study by Bal- 
bus & Hawley (1992a) suggesting that the Oort constant 
A is the maximum growth rate for such shearing instabil¬ 
ities. 

The instability still exists when the magnetic field is 
strong, but the growth rate decreases as This in¬ 

stability may play an important role, not only in weakly 
magnetized accretion discs, but also in the gaseous disc of 
galaxies, where the magnetic pressure is comparable to the 
thermal pressure: if a 1 and A/^l = —0.5, the e-folding 
time of the shearing instability is about one half of the 
rotation period for a to = 1 perturbation at i? ~ lOkpc. 

We can notice here an important difference with the 
axisymmetric instability of a disc embedded in vertical 
magnetic field: in this latter case, the instability criterion 
again compares the differential force to the magnetic ten¬ 
sion: AAQ. + klV'^ < 0. In particular, this condition implies 
that the instability cannot exist if the magnetic energy 
is comparable to or larger than the thermal energy. In¬ 
deed, the vertical wavelength of the perturbation must be 


shorter than the scaleheight of the disc, so that axisym¬ 
metric instability requires (Balbus & Hawley 1991): 

Va < —~ Max(as, Va). (33) 

TT 

By contrast, the azimuthal wavelength is only bounded 
by the circumference of the disc, so that the non- 
axisymmetric instability we are studying occurs as long 
as: 

Va <~ Max(as, Va)^- (34) 

TT H 

The shearing instability might therefore be efficient for 
magnetic fields exceeding the equipartition value, in the 
outer parts of the disc. This might be an important prop¬ 
erty of the instability, especially if it were related to a 
dynamo mechanism. 


4-2. Comparison of the growth rates 

It is interesting at this point to compare the characteristic 
length scales and growth rates of the optimal shearing and 
Parker instabilities. Both occur at large azimuthal scales, 
i.e. at low ky. The Parker instability requires ky < fcp at 
\kx\ oo (where fcp is defined in Eq. (|A8[)), whereas the 
shearing instability requires fcy < fc^ at [fc^j —> oo (where 
kA is defined in Eq. (p^). 

As an illustration, FigTO (top) shows the shearing instabil¬ 
ity occurring at an azimuthal wavelength which is stable 
against the Parker instability {kA > fcj/ > fcp). It was ob¬ 
tained by integrating numerically the differential system 
&■ The bottom graph shows the corresponding slow- 
wave solutions of the WKB dispersion relation. The two 
solutions (corresponding to slow waves propagating in op¬ 
posite directions) are not exactly symmetric with respect 
to the Lo axis, showing an effect of rotation which was 
already noted by Shu (1974). The two solutions merge 
at Kx — —kz where they become unstable. The numeri¬ 
cal solution shows that, as a result, when one starts with 
a “pure” solution on the top branch, one gets at posi¬ 
tive Kx a mixture of the two waves. This illustrates the 
linear coupling, which occurs because the WKB criterion 
is not strictly satisfied where the two slow MHD waves 
merge: the change of polarization and frequency between 
the oscillation phase and the growing phase is too fast (the 
shear, corresponding to a flat rotation curve, is strong: 
—A/Vt = .5). Figureillustrates this rapid change. 

We have plotted in Fig. ^ the range of magnetic pres¬ 
sures and rotation profiles for which, for any azimuthal 
wavenumber, the shearing growth rate is always higher 
than the Parker growth rate (contour levels > 1), and con¬ 
versely (contour levels < 1). For instance, along the con¬ 
tour level labeled “2”, the shearing growth rate is twice as 
large as the optimal Parker growth rate, at the wavenum¬ 
ber ky which is optimal for the Parker instability. 
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A “realistic” galactic disc (a ~ /3 ~ 1, —A/Q, ^ 

.5 to .75) would lie in the zone where the strengths of 
the two instabilities are comparable. On the basis of the 
comparison, in the linear regime, of the growth rates of 
these two instabilities, it is not possible to favour one par¬ 
ticular vertical wavenumber. This contrasts with the case 
of a uniformly rotating disc, where the maximization of 
the growth rate of the Parker instability over all possible 
wavenumbers would select the mode with kz = 0 (Parker 
1966, Zweibel & Kulsrud 1975). 

Non-linear effects will play a major role in selecting the 
vertical wavelength of the most unstable perturbations, 
but this is beyond the scope of the present paper. 

In an accretion disc with Keplerian rotation {—A/fl ~ 
.75), the growth rate of the shearing instability is higher 
than the Parker growth rate at any a < 2. In such a disc, 
the fastest growing perturbations, in the linear regime, 
have a short vertical wavelength and are essentially hori¬ 
zontal. 

We must bear in mind that this comparison is quan¬ 
titatively dependent on the relationship assumed between 
the scale height H and the rotation frequency fl (Eq. (^). 
Figure ^ is, however, qualitatively relevant. A lower scale 
height would slightly favour the Parker instability. On the 
other hand, if the scaleheight of the magnetic pressure 
were higher than the scale height of the gas, the Parker 
instability would be weakened, while the shearing insta¬ 
bility would remain unchanged. 


Fig. 3. With kA > ky > kp, the perturbation is stable 
against the Parker instability at \Kz;{t)\ > \kz\, and unsta¬ 
ble to the shearing instability at \Kx{t)\ < \kz\. The evolu¬ 
tion of the displacement vector (top figure), is in very good 
agreement with the calculation of the slow MHD root of 
the WKB-dispersion relation (bottom figure). Parameters are 
—A/Q = .5,kz = 10., a = 0.1,/3 = 0. The three components 
{x,y,z) of the displacement vector ^ are displayed. Wavenum¬ 
bers are in units of H~^, frequencies in units of VaIH. 


In the transition between these two regions, the choice of 
the azimuthal wavenumber determines which instability 
is faster. Introducing a cosmic-ray gas favours the Parker 
instability only: the pure-shearing instability does not de¬ 
pend on the cosmic-ray pressure, because the displace¬ 
ments involved are essentially horizontal and the differen¬ 
tial force does not act on the massless cosmic-ray gas (see 
Eq. (il)). 


Fig. 4. Comparison of the growth rates of the shearing and 
Parker instabilities, depending on the strengths of the magnetic 
field and differential rotation. The contour levels give a con¬ 
servative estimation of the ratio ujs/ujp\ it is minimized (resp. 
maximized) over every ky when ws/ujp > 1 (resp. ujs/wip < 1). 
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This comparison of the two isolated optimal instabili¬ 
ties is of course only a first estimate of their relative impor¬ 
tance. A refined study would require taking into account 
the temporal evolution of the radial wavenumber, and not¬ 
ing how the instabilities may be simultaneous or mutually 
exclusive. 

We already know that the Parker instability is favoured 
by large radial wavenumbers and vanishing vertical ones, 
whereas the shearing instability favours large vertical 
wavenumbers and limited radial ones. The next sections 
are devoted to a further study of this point. 

5. The Parker-Shearing instability 

5.1. The simplified dispersion relation for slow MHD 
waves 


This indicates that the dispersion relation ( |^ ) gives, for 
hnite values of Kx{t), the behaviour of the two slow MHD 
frequencies. 

Now we study the stabilizing or destabilizing influence 
of the differential force on the slow MHD branch. 

5.2. Stabilization of the Parker instability by differential 
rotation 

When k1 < 1/4, a strong differential force > ky) 
impedes the Parker instability at low \Kx{f) \ even though 
ky + kl < fcq. This stabilization was found in [FT] for kz = 
0, and is illustrated in Fig. || for kz = 0.25. Nevertheless, 
this stabilization lasts a few shear times only. 


As mentioned by Balbus & Hawley (1994), the shearing in¬ 
stability is the consequence of the differential force rather 
than of the shearing of the flow. In particular, this local 
instability would not exist in a shear flow between two 
parallel planes, as can be easily checked by canceling the 
rotation terms (Coriolis and differential forces) and keep¬ 
ing the shearing of waves (time-dependent wavenumber). 

Nevertheless, the sheared motion is responsible for the 
transience of the instability, as well as a linear coupling 
of the waves at low \Kx{t)\, as was shown in [FT]. We 
concentrate here on the nature of the Parker and shearing 
instabilities, and avoid the complications due to this cou¬ 
pling by considering the case of low shear {A/Vi 0). 

In order to study the effects of the differential force, 
we keep it comparable with the magnetic tension, i.e. 

klVl ~ 4AD. (35) 

This leads us to consider the case of a vanishing magnetic 
field (a —> 0) so as to keep: 

kA ~ ky, (36) 

~ 1/a » 1. (37) 


5.3. Destabilization of the slow MHD mode by differential 
rotation 


When kf > 1/4, the discriminant of Eq. ( p^ ) is negative 
if Kx{t)’^ < Kf, with: 


= {Tkl + 1 ) 


kl{kl-kl-kl) + k\{kl-l/A) 

Aklikl + 1 / 2 ) + kl- 1/2 


(41) 


The differential term k\ clearly contributes to destabilize 
the slow MHD mode, resulting in the shearing instability. 


When kz ky,Kx,l/2, we recover the pure shearing 
instability studied in Sect. 4: 

(42) 

with the usual criterion for instability: 

kA>ky^ 4AH -k k^Vl < 0. (43) 

For ky = k\/2 we recover the optimal growth rate: 

|ws| ~ \A\. (44) 


Hence we use the same WKB-approximation as in [FT], 
but include here a non-vanishing vertical wavenumber kz. 
Assuming that Cj = Oil/Tl), we can write the leading 
terms of the WKB-dispersion relation (Bll) and obtain: 

(kl + {Co^f + 2Kjiykz{n0j) - kyR = 0, (38) 

with: 


R = klikl - kl) + kim - kl) + kl/A - Kl{kl - kl).m 

When kA = 0, we recover the classical criterion for Parker 
instability at = 0: 

.2 


5 . 4 . Efficiency of the shearing instability at kz > 1/2 

Eq. (p^ is useful in determining the overall strength of the 
transient shearing amplification. When the discriminant of 
Eq. ( psi) is negative, the complex root is: 


~ kykzKx r., I . 

Hw = -— 1 ± i 

kl + l/M 


1 + fc + 


—V 

klKli 


(45) 


We define the complex amplification factor ip so that the 
amplitude of the wave after the transient amplification is 
multiplied by exp((p): 


<p= iu!{t)dt = 


Va 


r+Ks 


0 Ah TJ I 


ioj{Kx)dKx. 


,2 


,2 


(46) 
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5.5. Comment on the transient character of the Shearing 
instability 

We wish to emphasize that the shearing instability, al¬ 
though transient and occurring preferentially at low kx, is 
essentially local and does grow on a sufficiently short time 
scale that its transient nature poses no physical limit to 
its growth: this contrasts with e.g. spiral density waves, 
which travel in the disc and can experience amplification 
only for a limited time (as they are reflected from coro¬ 
tation), unless the boundary conditions at the disc center 
allow them to be reflected back toward corotation. 


Indeed the time required for the growing perturbation 
to travel, vertically or radially, along a distance compa¬ 
rable to the scaleheight of the disc is much longer than 
the typical e-folding time for realistic values of differential 
rotation. This can be checked by derivating the dispersion 
relation (B14) with respect to kx (or kz) to obtain the 
radial (or vertical) group velocity: 


Vgx _ 1 du! 

IT ~ Hdk^ 


n 

Hkz 


«l^l 


(49) 


Thus perturbations can be considered to grow where they 
were created. 


The transient character of any instability may be an 
obstacle to its development only if the instability does not 
have enough e-folding time to produce a significant effect: 
we have shown that the shearing instability does not suf¬ 
fer this limitation if we allow for a short enough vertical 
wavelength. 


Fig. 5. Transiently stabilized Parker instability. The shearing 
instability criterion is satisfied {ky < Ua), the differential ro¬ 
tation is stabilizing because < 1/2. Parameters are a = .4, 
/? = 0, —A/Q, — .1, ky = kz = 0.25. Consequently, fcq = 0.67, 
and kA ^ kv = 0.84. Here, the stabilization lasts about 4 shear 
times 

The imaginary part of (p vanishes because Im((I>) is odd in 
Kx, and we obtain: 


Incidentally, it has been shown in the context of non- 
magnetized discs that any local instability is stabilized by 
the effects of viscosity in a differentially rotating disc (see 
Korycansky (1992) or Dubrulle & Knobloch (1992)). This 
stabilization arises when the sheared radial wavenumber 
increases towards infinity, i.e. when the radial wavelength 
is smaller than the typical scale for dissipative processes. 
Nevertheless, such a mathematical limitation of the lin¬ 
ear regime has no physical implication if the “transient” 
instability has had enough e-folding time to produce per¬ 
turbations which lead to the non linear regime. 


mPci-Py-Pz) + ~k\{Pz-y^) 

kl (4fc2(fc2 + l/2) + fc2_l/2)l/2 

The amplification factor is an unbounded increasing func¬ 
tion of kz- For large values of kz, it scales as: 


,7r kz k\-kl 

2 (fc2 + i/2)i/2 k\ 


(48) 


The transient destabilization may be arbitrarily long if the 
vertical wavenumber kz is large enough. 


As an illustration, we know that the Parker instabil¬ 
ity is stabilized by ambipolar diffusion for scales < Ipc 
in the galaxy (see Cesarsky 1980): this leaves more than 
100 e—folding time for the Parker instability to occur on 
larger radial scales. 

Similarly, the shearing instability is limited by the min¬ 
imum vertical scale at which the ideal MHD equations are 
valid. It can be estimated, in a galactic disc, as the same 
critical scale used in the radial direction for the Parker 
instability. Although transient by definition, the shearing 
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instability is not physically more transient than the Parker 
instability. 


Fig. 6. Shearing instability at \Kx\ < Kg ~ 20, and Parker 
instability at \Kx\ > Kp ~ 40, separated by an oscillating 
phase. Here the shearing growth rate is lower than the Parker 
one. Parameters are a = 1, = 0, —A/Q, = .25, and ky = .57 

is optimal for the Parker instability kz = 20. We have then 
kA = kp = 1. and kq = .87 


6. Transition between the Parker and the shearing 
instability 

6.1. General evolution of the unstable slow MHD wave 

The dispersion relation demonstrates that both the 
shearing and the Parker instabilities belong to the same 
slow MHD branch, thus confirming the qualitative argu¬ 
ment of Sect. 3.3.: the Parker instability occurs when the 


Fig. 7. Shearing instability at \Kx\ < Kg ~ 25, and Parker 
instability at \Kx\ > Kp ~ 65, separated by an oscillating 
phase. Here the shearing growth rate is larger than the Parker 
one. Parameters are a = 1, /3 = 0, —A/Q. = .5, kz = 20, and 
ky = .89 is optimal for the shearing instability. Then kA = 1.41, 
kp = 1., and fcq = .87 

polarization of the slow MHD wave is vertical {K^ff) <C 
— \kz\), and is stabilized when the polarization rocks to¬ 
wards the plane of the disc {\Kx{t)\ < kz). The differen¬ 
tial force then acts, in the horizontal plane, on the radially 
displaced gas in the same destabilizing way as did the gra¬ 
dient of magnetic pressure, in the vertical plane, on the 
vertically displaced gas. While the shearing goes on, the 
polarization rocks back towards the vertical plane, lead¬ 
ing to the Parker instability again at K^it) ^ kz. This 
behaviour is illustrated in Figs.(5-8). 

The relative growth rates of the two successive insta¬ 
bilities (Parker/shearing) simply depend on the strengths 
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of the magnetic field and differential rotation, as studied in 
Sect. 4. The two instabilities can follow one another with¬ 
out any interruption of growth as Kx{t) varies (Fig. ^), or 
they can be separated by a phase where the slow MHD 
mode is oscillating as in Figs.(6-7). 


6.2. The case of weak differential rotation 


Even without differential rotation, the Parker instability is 
naturally stabilized at low kx when the vertical wavenum¬ 
ber is high. This appears on the classical criterion for the 
Parker instability kf + kf < 


kq at kx = 0. This defines 
a radial wavenumber Kp such that stabilization occurs at 
Kp. It can be calculated from the classical dis- 


\kx I ^ 

persion relation without differential rotation ( B12 ). An 
easy estimate of Kp is possible in the absence of rotation 
(O = = 0), for large values of k^: 


K 


kl 


1/2 


'P I ^2 1.2 

\ 


k,. 


(50) 


The shearing instability occurs at low \Kx{t) \ only, and 
is stabilized for \Kx{f)\ > Kg even when the Parker insta¬ 
bility is not present (see Balbus & Hawley 1992b). The 
value of Ks depends on the intensity of the differential 
rotation as in Eq. &■ Kg is negative if kA < ky, and 
it naturally increases when the shear parameter increases. 
We obtain, for kz 'i^ky-. 



A weak differential rotation is therefore likely to imply 
Kg < Kp. 


6.3. The case of realistic differential rotation 

We can derive a criterion for the presence of an oscillating 
phase of transition between the two instabilities, for large 
values of kz. According to our previous estimate, such a 
transition occurs at Kx kz. The WKB-approximation 
is then possible even if the differential rotation is strong 
(see Appendix B). 


If an oscillating transition occurs, the two slow MHD 
frequencies Cj{Kx) are purely real, and satisfy the following 
dispersion relation: 


Qa + 



2 

A — 21215(1 ey P)ky 



= 0 . 


(52) 


The two frequencies oj{Kx) pass through an extremum 
(see Fig. ^ and Fig. ^ at lo±{K±) such that: 


doj ^ ^ K± Qu}±{l + a + P)ky 

dkx kz A 


Fig. 8. Parker-Shearing instabilities occurring successively: we 
can clearly see the two instabilities occurring at the same time 
while Ks > \Kx{t)\ > Kp, with Kg ~ 50 and Kp ~ 20; this 
appears as a distinct bump in the instant growth rate (bot¬ 
tom figure). The shearing growth rate is slightly larger than 
the Parker one. Parameters are a = 1, 13 = 0, —A/Q = .75, 
kz = 20, and ky = .57 is optimal for the Parker instability. 
Then kA = 1.73, kp = 1., and fcq = .87 


In the limit ot kz 3> 1, the two extremal frequencies 
uj± are thus solutions of the polynomial of degree 4 in d)^: 

AQA = n^d}^{l + a + pfkl (54) 

The limiting case of a vanishing period of oscillations cor¬ 
responds to (5+ = d)_ and Ar+ = K-. This condition 
means that the polynomial (^) has a double root. It can 
be translated into a condition on its coefficients, depend¬ 
ing on ky, a, (3, A/Vl. It is solved numerically and displayed 
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on Fig. We can stress the following properties, which are 
independent of (a,/3): 


Fig. 9. Range of azimuthal wavenumbers leading to an os¬ 
cillating transition between the instabilities {Ks < Kp), or a 
phase of simultaneous instabilities {Ks > Kp), depending on 
the strengths of the magnetic field and differential rotation. 
The dotted line corresponds to Kg — Kp. In the top figure, 
/3 = 0 and a = 1. At —A/Q, = 0.5 (bottom figure), the two in¬ 
stabilities are never simultaneous. The two dotted curves cor¬ 
respond to the two critical wavenumbers ky leading, for any 
strength a of the magnetic field, to Ks = Kp. Any other ky 
implies Ks < Kp. 

(i) If the rotation curve decreases more slowly than 
the flat rotation curve {—A/Q. < .5), the two instabilities 
will always be separated by an oscillating transition {Kg < 
Kp). If —A/Q. is close enough to the critical value (.5), the 
duration of the transition can be arbitrarily short: there is 
a unique azimuthal wavenumber ky leading to Ks = Kp . 


Even in this case, the transition is slow because the growth 
rate still vanishes between the two instabilities (see Fig. 
as an illustration). 


Fig. 10. Optimal transition between the Parker and the shear¬ 
ing instabilities when —A/Q, < .5: the growth rate vanishes 
between the two instabilities. Here, a = 1, /3 = 0, —= -4, 
fez = 20, and we choose ky = .48 according to Fig. H so that 
Ks = Kp{~ 35). Then fcA = 1.26, fcp = 1., and fcq = .87 


(ii) If the rotation curve decreases faster than the flat 
rotation curve {—A/Vl > .5), the two instabilities can fol¬ 
low each other continuously, and even be simultaneous: 
Ks > Kp as in Fig. || where the rotation is Keplerian 
{-A/n = .75). 

(hi) If the rotation curve is exactly flat (—A/O = .5), 
Fig. ^ (bottom) shows that there are two critical values 
of ky allowing Ks = Kp, for any strength of the magnetic 
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field. This demonstrates that the lowest value of —A/Vl 
permitting Kg > Kp is precisely —A/Q, = .5, for any a. 

7. Discussion and conclusions 

Our study reveals that the WKB approximation is a very 
powerful tool to study analytically the effect of differential 
rotation on the slow MHD waves. Our numerical compu¬ 
tations are in very good agreement with the WKB disper¬ 
sion relation, even when the WKB criterion is not strictly 
satisfied. We summarize here our main conclusions: 

(i) The shearing instability comes from the non- 
axisymmetric destabilization of the slow MHD wave by 
the differential force. However, a parallel can be drawn 
with the axisymmetric shearing instability of the Alfven 
mode in a disc embedded in a vertical magnetic field. 

(ii) The Parker instability and the shearing instability 
both correspond to the destabilization of the same slow 
MHD mode. The occurrence of each instability can be 
understood by looking at the natural polarization of the 
slow MHD mode, depending on whether their radial wave¬ 
length is longer or shorter than their vertical wavelength. 

(iii) Although transient, the shearing instability may 
last an arbitrarily long time if the vertical wavelength is 
short enough. In this respect it is a true local instabil¬ 
ity. Moreover, it still exists for a magnetic field strength 
exceeding the equipartition value. 

(iv) The action of differential rotation is opposite de¬ 
pending on whether the vertical wavelength of the slow 
MHD wave is longer or shorter than the scale height of 
the disc. Perturbations with a long vertical wavelength 
{kz < \/2H) are transiently stable against the Parker in¬ 
stability, during a few shear times. On the contrary, slow 
MHD waves with a short vertical wavelength are destabi¬ 
lized by differential rotation almost as long as their polar¬ 
ization is in the plane of the disc, i.e. while \Kx{t)\ < kz- 
The Parker instability naturally occurs when the po¬ 
larization of the slow MHD wave is vertical, i.e. when 

\Kz,{t)\ > kz- 

(v) A refined look at the transition between the Parker 
and shearing instabilities shows that the intensity of the 
differential rotation is the crucial parameter. If the rota¬ 
tion curve decreases more slowly than the flat rotation 
curve {—A/Q. < 1/2), the growth rate will always van¬ 
ish between the two instabilities. If the rotation curve de¬ 
creases faster than the flat rotation curve {—A/Q. > 1/2), 
the two instabilities can occur simultaneously: in discs 
with a strong differential rotation, the most unstable per¬ 
turbations may not be the ones with kz = 0, because of 
this shearing instability. 

Astrophysical and observational consequences of these 
results are beyond the scope of this paper, which leaves 
many questions open: observing shearing instabilities 
might prove even more difficult than looking for the sig¬ 


nature of Parker instabilities, because here one is dealing 
with small vertical wavelengths, on which the gas disc and 
the field are hardly homogeneous. This might actually be 
seen as an obstacle to the development of the instability, 
or as a consequence of it (see the discussion in Zweibel & 
Kulsrud 1975, Parker 1975b). 

Furthermore, since we are limited to linearized analysis, 
we cannot answer the question of the ultimate fate of 
these waves. If they reach very large amplitudes before the 
shear takes them to very large they might more easily, 
through the non-linear effect of self-gravity on the result¬ 
ing clumps of gas, lead to Jeans collapse of these clumps. 
On the other hand if they reach large their ultimate 
fate would be more likely a cascade towards small-scale 
turbulence. 

One should also remember, in such a discussion, that in 
galaxies the gas disc is swept and shocked by the spiral 
arms on a time scale of the order of a rotation time, i.e. at 
best a few amplification times, so that the coherent growth 
of a Parker or Shearing instability on a longer time scale is 
rather dubious. The conclusions would obviously depend 
on assumptions on the initial state, i.e. the amplitude of 
fluctuations generated in the disc and ready to be ampli¬ 
fied by these instability mechanisms. 

On the other hand self-gravity might play a quite dif¬ 
ferent role, which we consider as most promising for gen¬ 
erating strong Parker-like perturbations in the disc: here 
and in [FT] we have shown that differential rotation re¬ 
sults in a linear coupling of waves, i.e. a wave on a given 
branch of the dispersion relation evolves into a mixture of 
the three MHD waves (slow and fast magnetosonic, and 
Alfven waves). If we had included self-gravity in this dis¬ 
cussion, the fast magnetosonic wave would have appeared 
as the classical spiral density waves of galactic discs, modi¬ 
fied by inclusion of the magnetic field (as was done, in the 
case of a vertical field, by Tagger et al. 1990). But one 
can conclude that a spiral wave does, by this mechanism 
of linear coupling, generate Alfven and slow magnetosonic 
waves, i.e. Parker rather than shearing instabilities since 
this occurs at kz = 0. We know that the spiral wave is 
strong (the density contrast is of the order of 100% in 
the gas), and the coupling is also strong since the rota¬ 
tion curve has —A/H ~ .5. This means that the resulting 
Parker perturbation will also have a very large amplitude, 
even before it has started to grow by Parker’s mechanism. 
We consider this as the most promising way of generating 
strong Parker perturbations in the disc, and will discuss 
it in future work. 
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A. Appendix: Laplace Transform of the solution 
with shear 


The polynomial A, already introduced in [FT], identifies 
with the asymptotic dispersion relation at kx —^ oo: 


We define the dimensionless rotation parameter and 


epicyclic frequency as: 



(Al) 

P = 

Vl 

P = A^. 

(A2) 


In order to reduce the degree of the differential system 
(pd]), we introduce the convenient function: 




1 + 20 
2 a 



+ 


1 + 2a 


+ 



1 + 2/3 \ 

2 + 4ay 



(A3) 


The only time-dependent parameter, Kx, now multiplies 
the functions ^x and only. Thus we can write two equa¬ 
tions without any A^j-multiplication, and two other equa¬ 
tions expressing Kxix and Kx^w in terms of ^w,^x,^y,iz- 
These latter are written: 


k^x^x 


— — + P - P 
vl dt^ ^ y ^ 


A 
-2a 
l-b2a 


_ H d - 


(A4) 




l-b2a 


Cj/ ( ^2 + 


, 1 -|- 2/3 , / A r\ 

■ *- - ^2(A5) 

2 -b 4a ' ^ ^ ’ 


Performing at this point the Laplace transform 
changes the A 2 ,-multiplications into w-derivatives, and 
time-derivatives change to simple ^-multiplications, with 
some additional w-polynomials whose coefficients de¬ 
pend linearly on the initial conditions ^{kx,t = 0) and 

§l£.{kx,t = 0 ). 

In the following, we use the dimensionless variable Cj = 
Hco/Va- 

The first couple of equations, independent of Kx, relate 
i^y,^z) to {^w,^x) in a simple algebraic manner: 


A 


— 

f Qa 

Qb \ 

6 

\qc 

Qd ) 


Qy 

Qz 


(A6) 


A = (1 -b 2a)u;^ - 2 |(1 -b a)ky -b (1 -b a -b /3)^^^— 

+ ~kl{Py-kl.)- (A7) 

The maximum azimuthal wavenumber kp allowing the 
Parker instability at kx oo is defined as: 


kp = {1 + a + (3) (A8) 

Equations (A4)-(A5) may be expressed as a differential 
system of second order on the functions {^w,Cx), using 

Eq. (@. 


tt duj 


II 

IvSl 

fPa 

Pb \ 


\Pc 

Pd ) 




(A9) 


where {pa,Pb,Pc,Pd) are simple polynomials of A whose 
coefficients are independent of kx- 

Pyj and Px are two polynomials of A of respective degree 5 
and 3, whose coefficients depend linearly on the six initial 
conditions 



^ ^ ^ 
dP""' dpy' dP" 


{kx,t = 0 ), 


(AlO) 


independently from the kx and parameters. 

The coefficients of the matrix in Eq. (|A^) are given by: 


Pa = V,kyUj 


Ay -b (1 -b a -b /3) 




Pb = -i{{P - Py 
Pc = -iQ2{P) - 

Pd = —tlkyUJ < P 


- P)AiP) - PPyiP -Py+kl)} 
~ (1 + 2 a)a;^}fc 2 

-kl + {l + a + f3) - ikx'^ I 


The determinant of (^) is PaPd —PbPc = A{Q^ — k^QA)- 
Defining Aq = fcp — 1/4 leads us to write the polynomials 
(52, Qa and Qa as: 


Qy and are polynomials of Cj whose coefficients de¬ 
pend linearly on the initial conditions. The polynomials 
[qa, Qb, Qc, qd) are independent of kx, and constitute a ma¬ 
trix which is invertible as long as A yf 0 and A yf 0: 


Q_a - Ay 

qb = /Dw 


— ky + [1 + a + P) (^kz -b —^ , 

|(1 -b 2a)(Ay — P) -b (1 -b a -b 0) 


Qc — 2^ky -b (1 -b 2P)P} + kz{ky — (1 -b 2a)A^}, 
qd = -flPy{l + a + 0). 


Q 2 = 2aP - (1 -b 2a) ^Ay -b -b ky[ky - Aq)(A11) 

Qa = (P -Py- P)Q2iP) - PyP[P PpQ- Py){kl2) 

Qa = (1 T 2(x)P — ti;^{2(l -b ct)ky -b (1 -b 2a)At^} 


+Ay(Ay — Aa)- 


(A13) 


The differential system ( [A^ ) can easily be transformed, by 
substitution, into the second order differential equation: 


f d'^ dlogfip d 
\ dP doj dtu 


+9[pH[kx,p = h[kx,p, (A14) 






















16 


T. Foglizzo et al.: The Parker-Shearing instability in azimuthaly magnetized discs 


where h depends linearly on the initial conditions, and 
(/,(?, h) are singular where A(ci)^) = 0. 

We derive for 


dlogfx 91ogA/pc ^ Pa+Pd 


doj 


duj 


k\ky A 


I^A^y 


n A du Pc 


W{iv^)^g{co^)- 


152 log/((D^) l/aiog/(A2) 




4 \ duj 


1 9^ log W fid log W 


4 9a;2 


4 doj 


< 9, 


^±(A) = 


1 


ipwy 


■ exp 


±i j W^/2(d;'^)dA'j . 


= - 


Vt 


k\ky 


The solution of the differential equation (A14) is then 
written as: 


1 cip+oa 

i{kx,t) = — / (/r+e*'"+ - 

Jip—oo 


W{P) = 




1 . 4 : 1.2 
t^A^y 


Kx{t) - 


Pa +Pd 

2iA 


The function WfjP) defined in 

^2 r 


w = - 


(A15) 


fc2fc4A2 


PcPb 


Pa - Pd 


|) can be written: 

1 1 92 A 

f 2 ^ 


^Ipc 


\Q 3 -klQA k\kyPc d Pd\ 

) = TTTF i -X- ) ■ (A16) 


d Pa-Pd ^ f d , A^ 

-ils5“>«yj 


B. Appendix: WKB approximation 

We concentrate here on the function ^x, without specifying 
the subscript x for the sake of clarity. The homogeneous 
differential equation associated to Eq. ( A14| ) may be writ¬ 
ten in its canonical form: 

+ = (B1) 

where W(d)2) is a fractional function of di2 defined as: 


2kyk\A duj Pc 

The general dispersion relation becomes: 

kyk\Apc d Pa-Pd 


{pa-iKxA){pd-iKxA) -pcPb 


PMA2 r g2 


A 


2 C 1 duj Pc 


d 


We recognize the first two terms as the determinant of the 


.(B2) 


We have shown in [FT] that if the WKB criterion is sat¬ 
isfied, 


system (A9). Some additional terms must be taken into 
account, especially in the limit of both vanishing shear 
and magnetic held (Sect. 5). 

It can be divided by A to obtain: 

Qd. — QaP; — Akl + 2fiA(l -I- a -I- f3)kykxkz 

+kyk\ i"3(1)2 _ ^2 ^ /3) — 


(B3) 


-2kyk\uj‘^ (1)2 - /c2 _|_ (1 _|_ Q, 


d log Pc 

2 a 


the solutions of the differential equation (^) can be ap¬ 
proximated by a linear combination of the two indepen¬ 
dent functions: 


klk\A 52 , A / 5 , A 


Pc 


dCo^ 


= 0 . 


(Bll) 


(B4) 


According to the dehnition (A15) of the function /, we 
can write: 

Let us dehne the phase lE'±((I>,t) in a dimensionless form 
as: 


Both the WKB criterion (^) and the dispersion relation 
(Bll) become much simpler in the two cases we investi¬ 
gate in Sect. (4-6): 

(i) If the shear vanishes totally, we will obtain the usual 
dispersion relation: 

Qs - QaP, - ^kl + 2flCj{l -b a -b d)kykxP = 0. (B12) 


(ii) If, as in Sect. 4 and Sect. 6, fcz 
obtain: 


9{P) 


L!2 




klk\A 


QA{P)k 


00 , we simply 


(B13) 


W{P) - 


L !2 


fc2fc4A2 


(-■ipb)(A:2 - (1-b 2a)(D2)fc2. (B14) 


Thus the criterion (B3) is satished in the limit of large 


(B7) 


where the two functions p±{Cj) vary slowly if the WKB 
criterion is satished. The dispersion relation is therefore: 


(B 8 ) 


vertical wavenumbers. 

In Sect. 4, the radial wavenumber is bounded {\Kx{t)\ <C 
jfczj), and the dispersion relation (Bll) becomes Qa{P) = 

0 . 

In Sect. 6 , \Kx{t)\ ~ \kz\ ^ 1/H, and the dispersion rela¬ 
tion (Bll) becomes: 

QT+(i^) A-2flA{l+a+l3)~ky = O (B15) 
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(iii) If, as in Sect. 5, a —> 0 and ~ cte, then: 

Wiu^) ~ 5(^2) ~ 0(0^)^ (B16) 

and the WKB criterion is again satisfied. The simple dis¬ 
persion relation (^) describing the slow MHD frequencies 
for finite values of kz,Kj. comes from the leading terms, 
when wfl ^ 1, of: 

Qs — Q-\- oi -\- P^kykxkz 

-\-kyk\ — ky {1 Oi ^ = 0- 


C. Appendix: polarization of a slow MHD wave 


Here we recall the {kx, )-dependence of the polarization 
of a slow MHD wave in absence of shear. Strictly speaking, 
the parameter kA ought to be zero in what follows. 

If to is determined by the dispersion relation (B12), the 
equations (|A^ and (|A6D will imply: 


kzQA - H(h(I -I- a -I- P)kykx + iRi 


(Cl) 


fca;A — Hd;(l -I- a -I- I3)kykz + ii?2 
where are independent of {kx,kz)'- 

Q {le - + (3klQ^ + 


i?2(w) = Hdifcj, 


u)^ - kl + {I + a + (3)-^ 


(C3) 


Let us study the polarization of the Alfven mode: we 
have seen in Sect. 4 that its eigenfrequency is a root of 
= 0 when the vertical wavenumber kz ^ oo, and 
we know from [FT] that it is a root of A(a)^) = 0 when 
kx oo. Unlike the magnetosonic eigenfrequency (fast 
MHD wave), it remains finite in both cases. Thus this cal¬ 
culation is valid for both the Alfven and the slow MHD 

mode. _ 

In the limit oi kx oo, the dispersion relation ( |B12 ) 
shows that: 


kxA ~ 2H(I!a( 1 -I- a -|- ( 3 )kykz. 

The polarization dH) is then: 

Cz — Hcja)! + Q: +/3)fci/ 

Cx HAa)! o: (})kykz *I?2(Aa) 


(C4) 


(C5) 


corresponding to a vertical mode. 

In the limit of ^ oo, the dispersion relation (HI) 
implies: 


kzQA ~ 2Ha;s(l -I- a -|- P)kykx- 
The polarization dH) becomes: 


^X 


+ a + P)kykx + iRi{(^s) 1 
—i7u^s(l a P')ky k^ 


(C6) 


(C7) 


and corresponds to a horizontal mode. 

If both kx,kz 3> fcy, the polarization (Cl) is simply: 



(C8) 


Thus the plane of polarization of the slow MHD (and 
Alfven) mode is naturally vertical or horizontal, depend¬ 
ing on the ratio kx/kz- 

These properties are specific of the slow magnetosonic 
mode in a stratified rotating disc. In a uniform plasma 
{H —> oo and H —> 0), the polarization of the slow magne¬ 
tosonic mode is indeed completely different, and is exactly: 




kx 


(C9) 


It can easily be shown that: 

(i) the property (C5) is an effect of stratification, and 
is valid only if kx Hk'^n 

(ii) the property ( P^ is an effect of rotation, and is 
valid only \i kz kyol/Vl. 
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